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ABSTRACT 

" Aims. Our aim is to determine the presence and location of the emission from polycyclic aromatic hydrocarbons (PAHs) 
I— I towards low and intermediate mass young stars with disks using large aperture telescopes. 

4^ Methods. VLT-VISIR N-band spectra and VLTTSAAC and VLT-NACO L-band spectra of 29 sources are presented, spectrally 
resolving the 3.3, 8.6, 11.2, and 12.6 /xm PAH features. Spatial-extent profiles of the features and the continuum emission 
^ have been derived and used to associate the PAH emission with the disks. The results are discussed in the context of recent 
$^ PAH-emission disk models. 

Results. The 3.3, 8.6, and 11.2 jj,m PAH features are detected toward a small fraction of the T Tauri stars, with typical 
Ci upper limits between 1 xl0~^^ and 5 x lO"^'^ W m~^. All 11.2 /xm detections from a previous Spitzer survey are confirmed 
with (tentative) 3.3 /um detections, and both the 8.6 and the 11.2 /um features are detected in all PAH sources. For 6 
I detections, the spatial extent of the PAH features is confined to scales typically smaller than 0.12-0.34", consistent with 
^ the radii of 12-60 AU disks at their distances (typically 150 pc). For 3 additional sources, WL 16, HD 100546, and TY 
CrA, one or more of the PAH features are more extended than the hot dust continuum of the disk, whereas for Oph 
IRS 48, the size of the resolved PAH emission is confirmed as smaller than for the large grains. For HD 100546, the 3.3 
OO A*ni emission is confined to a small radial extent of 12±3 AU, most likely associated with the outer rim of the gap in this 
O^l disk. Gaps with radii out to 10-30 AU may also affect the observed PAH extent for other sources. For both Herbig Ae and 
^ T Tauri stars, the small measured extents of the 8.6 and 11.2 fjxn features are consistent with larger (> 100 carbon atoms) PAHs. 

Key words. Stars: pre-main sequence - planetary systems: protoplanetary disks - Circumstellar matter - Astrochemistry 



1. Introduction tion of the T Tauri stars, 8% (Gccrs ct al. 2006), show clear 
^ polycyclic aromatic hydrocarbon (PAH) features, which is 
Mid-infrared (IR) spectroscopy from the ground, with the compared to the 54% detected in Herbig Ae/Be stars 
Infrared Space Observatory (ISO) and recently with the (^cke & van den Anckcr 2004; Habart et al. 2004). Geers 
Spitzer Space Telescope has revealed that low (T Tauri) (2OO6) argued that their low PAH detection rate is 
and intermediate mass (Herbig Ae/Be) pre-main-sequenc;e consistent with a 10-lOOx lower PAH abundance in the 
stars often show silicate bands in emission (e.g., Meeus ^-^^^ compared to the ISM. One difficulty has been to 
et al. 2001; Honda et al. 2003; van Boekel et al. 2004; p^perly identify PAH features in mid-IR spectra in the 
Kessler-Silacci et al. 2006) . Nevertheless only a small frac- presence of strong silicate features. The 11.2 /zm band can 
be confused with the 11.2 /im crystalline forsterite feature 



Send offprint requests to: V. C. Geers, and can also be blended with the broad amorphous sili- 

e-mail: vcgeer sOstrw . leidenuniv . nl cate feature whose strength and spectral width varies with 

* Based on observations obtained at the European Southern grain size. 
Observatory, Paranal, Chile, within the observing pro- 
grams 164.1-0605 (ISAAC May 2002), 074.C-0413 (NACO, The 3.3 ^im feature, obtainable through ground-based 

March/April 2005), 075.C-0420 (ISAAC August 2005), 077.C- studies, is another diagnostic of the presence of PAHs in 

0668 (VISIR/ISAAC April/May 2006) disks. This broad feature is comparatively isolated and is 
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expected to be tightly correlated with the 11.2 /zm fea- 
ture since both involve C-H vibrations, as confirmed by 
ISO data (e.g. Peeters et al. 2004). Confirming the pres- 
ence of PAHs in the few T Tauri sources where they have 
tentatively been seen is the first aim of this paper. 

The second aim is to determine whether the PAH emis- 
sion comes from the disk or rather from an extended rem- 
nant envelope of dust around the star. With Spitzcr spec- 
troscopy, the spatial resolution of ~ 3" at 10 /zm in the 
high-resolution mode corresponds to spatial scales of 300- 
900 AU for typical nearby star- forming regions {d = 100 
300 pc). Ground-based observations of the 3.3, as well as 
the 8.6 and 11.2 /xm bands with 8m class telescopes, pro- 
vide a spatial resolution of < 0.3" (at 10 /im), which is 
an order of magnitude higher than that of Spitzer. Recent 
ground-based, spatially resolved observations of a num- 
ber of Hcrbig Ac stars at typical distances of 100 150 pc 
have shown that the measured PAH features come from 
regions with sizes typical of a circumstellar disk (radius 
< 12 AU at 3.3 /7,m, <100 AU at 11.2 /im) (Habart ct al. 
2004; Geers et al. 2005; Lagage et al. 2006; Doucet et al. 
2007). Spatially resolved, mid-infrared spectroscopy was 
presented by van Bockcl ct al. (2004) for 3 bright Hcrbig 
Ae stars. They found that the 8.6, 11.2, and 12.7 /xm fea- 
tures are extended with respect to the continuum emis- 
sion, on scales of (several) 100 AU. Habart et al. (2006) 
presented VLT-NACO observations resolving the 3.3 /xm 
PAH feature above the continuum emission for 4 Herbig 
Ac/Be sources, with typically 50% of the intensity com- 
ing from within radii smaller than 30 AU. In contrast, 
observations of WL 16 by Ressler & Barsony (2003) show 
spatially resolved PAH emission with an extent of 880 x 
400 AU at 1% of the peak level. 

The above-mentioned observations support the recent 
modeling results by Habart ct al. (2004), Geers et al. 
(2006), and Visser et al. (2007), which indicate that most 
(~ 80%) of the spatially extended PAH emission comes 
from within a radius ^120-170 AU, with the exception 
of the 3.3 /xm feature for which they predict that half of 
the emission should originate from <50-160 AU from a 
typical Herbig Ae/Bc star, which depends strongly on the 
size of the PAH molecules. No such data exist yet for T 
Tauri stars, except for the unusual target IRS 48, a MO 
star for which Geers et al. (2007) measured PAH features 
with a radial spatial extent of ~ 75-90 AU. 

Determining the presence and location of PAHs in the 
T Tauri disks is significant for several reasons. Due to their 
high opacity at UV wavelengths, PAHs can be used as a 
tracer of the strength of the UV radiation field, while at 
the same time their presence in the inner disk can have a 
strong infiuence on the amount of UV that is received by 
the outer disk. Also, the strength of PAH features from the 
outer disk provides a tool to determine if the disk is fiat 
or fiaring (Dullemond & Dominik 2004; Acke & van den 
Ancker 2004; Habart et al. 2004). PAHs in the surface layer 
of the (outer) disk are an important heating mechanism 
of the gas through photo-ionization producing energetic 
electrons (Jonkheid et al. 2004). This in turn influences 



the outer disk chemistry and line emission. PAHs and very 

small grains (VSGs) also become an important site for H2 
formation when the classical grains have grown to jjm or 
mm sizes (Jonkheid et al. 2006). 

High spatial-resolution images taken in narrow band 
filters centered at particular emission features are one 
method of constraining the spatial extent of the emitting 
species, but a good discrimination between the contribu- 
tion of silicates and PAHs to the mid-infrared excess emis- 
sion is essential here. Long-slit infrared spectrometers such 
as ISAAC, VISIR, and NACO, installed on the 8m class 
VLT telescopes, allow for both spectrally and spatially 
resolved observations of disks around young (pre-) main 
sequence stars. 

In this article we present the results of two such 
surveys, carried out with the VLT-ISAAC and VLT- 
VISIR instruments. In addition, VLT-NACO L-band spec- 
troscopy is obtained on two sources, making use of adap- 
tive optics to obtain even higher spatial resolution. The 
aim of this work is to obtain a limited survey of the 3.3 
and 11.2 /xm features, with a focus on low-mass T Tauri 
stars with disks, to study the usefulness of the 3.3 /xm fea- 
ture as a PAH tracer compared with the 11.2 /im feature 
and to study the spatial extent of the PAH emission in 
the context of the disk. 

2. Observations and data reduction 

2.1. Source selection 

The sample contains 19 T Tauri stars and 10 Herbig 

Ae/Bc stars and was selected as follows: 17 T Tauri stars 
in the nearby star-forming regions Chamaeleon, Lupus, 
Ophiuchus, and Serpens were chosen from the sample ob- 
served with Spitzer Infrared Spectrograph (IRS) in the 
context of the "Cores to Disks" (c2d) Legacy program 
(Evans et al. 2003). All their sources with definite and 
tentative PAH detections (Geers et al. 2006) were cho- 
sen, including 4 Herbig Ae stars. In addition, as part of 
a backup program during pointing limited nights with 
strong winds, 4 additional young stars in the Serpens and 
Corona Australis star-forming regions, found in IRAC + 
MIPS imaging surveys, were observed with ISAAC. One 
source, HD 100546, was observed in an earlier ISAAC pro- 
gram. Three sources, IRS 48, WL 16, and EC 82, were 
added to the VISIR N-band program after serendipitous 
discovery of the 3.3 /im PAH band in the L-band spectra 
in an ISAAC survey of embedded sources (Pontoppidan 
et al. 2003). The flnal source list is given in Table 1. 

2.2. ISAAC L-band spectroscopy 

L-band spectroscopy was obtained with ISAAC, the 
Infrared Spectrometer And Array Camera, installed at 
the VLT Antu at ESO's Paranal Observatory in ChUe, 
in the nights of June 16, 2000 (HD 100546), August 9-14, 
2005, and April 17-18, 2006 in the low-resolution (A/AA 
= 600) spectroscopic mode in the spectral domain 2.8- 
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4.2 jj,m using a 0.6" x 120" slit. The telescope was oper- 
ated using a chop throw of 20" and a nod throw of 20". 
Because of chopping, these observations are only sensitive 
to spatially extended emission of at most 10". For flirx 
caHbration and telluric line correction, several standard 
stars were observed. The data were reduced using our own 
IDL routines, first described in Pontoppidan et al. (2003). 
The individual frames were corrected for the nonlinear- 
ity of the detector array and distortion corrected using 
a startrace map; bad pixels and cosmic ray hits were re- 
moved before co-adding, using a shift-and-add procedure 
to correct for telescope jitter. From the combined frames, 
both the positive and the two negative spectral traces 
were extracted and co-added. For correction of the tel- 
luric features, the extracted source spectrum was divided 
by the similarly extracted, standard star spectrum using 
an optimal small wavelength shift and an exponential air- 
mass correction between the source and the standard, thus 
requiring that the pixel-to-pixel noise on the continuum 
of the final science spectrum is minimized. Before this 
exponential airmass correction was applied, the detector 
and filter response curves were removed to obtain a spec- 
trum of the true atmospheric absorption. These curves 
were obtained by fitting an envelope to both the stan- 
dard star spectrum and a spectrum of the approximate 
atmospheric transmission and by taking the ratio of the 
two envelopes. Flux calibration was performed by divid- 
ing the science spectrum by an observed standard star 
and multiplying by a blackbody with the effective tem- 
perature of the standard star. Airmass correction was ap- 
plied by multiplying the flux of the science target by the 
factor io-0-4ExtL(AM.,t-AM3e) ^iicre we use for the L-band 
atmospheric extinction, ExtL = 0.08 magn. x airmass"^, 
the value provided by ESO on the ISAAC webpage^. The 
flux calibration is estimated to be accurate to 30%. The 
spectrum is wavelength-calibrated relative to the atmo- 
spheric transmission spectrum and is accurate to ~0.003 
fiia. Note that the observation campaign in August 2005 
suffered from poor weather conditions, with strong winds 
and variable seeing. Part of the nights were pointing lim- 
ited to the North due to strong winds, which led to a 
relatively large fraction of the backup program sample of 
Serpens sources being observed. A summary of the pre- 
sented observations is given in Table 1. 

2.3. NACO L-band spectroscopy 

L-band spectroscopy was obtained with NAOS-CONICA 
installed at the VLT-Yepun at ESO's Paranal Observatory 
in Chile, in service mode, in the nights of 2-5 26 March 
and 10 April of 2005 in the medium-resolution {R = 700) 
spectroscopic mode in the spectral domain 3.20-3.76 /xm 
using a 0.172"x28" slit. The data were reduced using our 
own IDL routines, following a similar procedure as for 
the ISAAC data. The spectrum was wavelength-calibrated 

^ http://www.eso.org/instruments/isaac/ 
imaging_standards.html 



relative to the atmospheric transmission spectrum and is 
accurate to ^0.008 /im. Flux calibration was performed 
using ISAAC spectro-photometry, estimated to be accu- 
rate to 30%. One source, WL 16, was observed with the 
long- slit of the spectrograph aligned in 2 settings, both 
perpendicular and parallel to the semi-major axis of the 
disk, which is at a position angle of 60 ± 2° (Ressler & 
Barsony 2003). A summary of the observations is given in 
Table 1. 

2.4. VISIR N-band spectroscopy 

N-band spectroscopy was obtained with VISIR installed at 
the VLT-Melipal at ESO's Paranal Observatory in Chile, 
in the nights of 3-7 May of 2006 in the low-resolution {R ~ 
350) spectroscopic mode in the spectral domain 7.7-12.5 
yum using a 0.75" x 32.3" slit. All sources were observed 
with the 8.5 and 11.4 /im settings to cover the main 8.6 
and 11.2 /xm PAH features. In addition, a few sources were 
observed with the 8.1 and/or 12.2 /xm settings to cover 
the red wing of the 7.7 /im, as well as the 12.7 /im PAH 
feature. Most of the nights were characterized by strong 
winds, pointing limited to the south and highly variable 
seeing. The telescope was operated using chop-throws and 
nod-throws of 8". Because of chopping, these observations 
are only sensitive to spatially extended emission of at most 
4". Several standard stars were observed for flux calibra- 
tion and telluric-line correction. The data were reduced 
using our own IDL routines, following a similar procedure 
as for the ISAAC data. Flux calibration was performed us- 
ing Spitzer IRAC 8 /xm photometry where possible. The 
uncertainty in the flux calibration was estimated to be 
30%. The spectrum was wavelength-calibrated using the 
ESO VISIR pipeline (v. 1.3.7) by comparison with atmo- 
spheric lines and is accurate to ^ 10~* fim. A summary 
of the observations is given in Table 1. 

2.5. Measuring spatial extent 

The FWHM of the spatial profile of the ISAAC spec- 
tra was derived from the co-added 2D spectral images, 
by fitting a Gaussian profile for each wavelength bin. For 
the ISAAC data, distortion correction was not performed 
because this correction was found to introduce a semi- 
sinusoidal pattern. The spatial extent of the PAH feature 
at 3.3 /xm was measured with respect to the extent of the 
disk continuum emission at 3.3 /zm, interpolated between 
the continua points at 3.1 and 3.5 /im. The extent in arc- 
seconds was converted to a radial extent from the center 
in AU for sources with known distances. 

Because of the variable seeing during a large part of our 
program, the atmospheric seeing was assumed to dominate 
the FWHM of the ISAAC and VISIR observations and 
the spatial blurring introduced by instrumental optics ig- 
nored. The standard star was assumed to be an unresolved 
point source and, where possible, its measured FWHM 
was used as a measure of the point-spread-function. Only 
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Table 1. Summary of observations. 



Target 


RA (2000) 


Dec (2000) 


Date 


Sp. Typo 


Dist. (pc) 


Ref. 


ISAAC L-band 


SX Cha 


10 55 59.73 


— 77 24 39.9 


17-04-2006 


MO 


178 


K06; W97 


SY Cha 


10 56 30.45 


— 77 11 39.3 


18-04-2006 


MO. 5 


178 


K06; W97 


WX Cha 


11 09 58.74 


—77 37 08.9 


19_ 04-2006 


K7-M0 


178 


K06; W97 


HD 98922 


11 22 31.67 


—53 22 11.4 


19-04-2006 


B9 


>540 


A98; H78 


HD 101412 


11 39 44.46 


—60 10 27.9 


19-04-2006 


B9.5 


160 


T94: A05 


HD 100546 


11 33 25.44 


— 70 11 41.2 


16-06-2000 


B9 


103 


A98: A98 


T Cha 


11 57 13.49 


-79 21 31.4 


18-04-2006 


G8 


66 


A93i A98 


IRAS 12635-7623 


12 57 11.73 


-76 40 11.1 


19-04-2006 


MO 




K06; - 


HT Lup 


15 45 12.86 


-34 17 30.6 


09-08-2005 


K2 


140 


H94; H93 


GW Lup 


15 46 44.73 


-34 30 35.5 


14-08-2005 


M2-M4 


140 


H94; H93 


SZ 73 


15 47 56.94 


-35 14 34.7 


13-08-2005 


MO 


140 


H94; H93 


GQ Lup 


15 49 12.10 


-35 39 05.1 


10-08-2005 


K7-M0 


140 


H94; H93 


HD 141569 


15 49 57.75 


-03 55 16.4 


13-08-2005 


B9.5-A0 


100 


AU04; Ai;04 


IM Lup 


15 56 09.22 


-37 56 05.8 


14-08-2005 


MO 


150-360 


H94: H93,K01 


RU Lup 


15 56 42.30 


-37 49 15.4 


10-08-2005 


K7-M0 


150-360 


H94; H93,K01 


DoAr 24E 


16 26 23.26 


-24 20 59.8 


16-08-2006 


KO 


125 


K06; G89 


Em* SR 21A 


16 27 10.28 


-24 19 12.7 


15-08-2005 


G2.5 


125 


P03; G89 


Em* SR 9 


16 27 40.29 


-24 22 04.0 


10-08-2005 


K5-M2 


126 


L99; G89 


Haro 1-17 


16 32 21.93 


-24 42 14.8 


14-08-2005 


M2.5 


126 


A93; G89 


V1121 Oph 


16 49 15.30 


-14 22 08.7 


10-08-2005 


K5 




VOO; - 


Wa Oph 6 


16 48 45.62 


-14 16 36.0 


11-08-2006 


K 




G07; - 


VV Ser 


18 28 47.86 


-1-00 08 39.8 


11-08-2006 


AOV 


259 


MOl; S96 


CoKu Ser G6 


18 29 01.23 


+00 29 33.0 


11-08-2005 


K3 


259 


C79; S96 


HD 176386 


19 01 38.89 


-36 53 27.0 


14-08-2005 


B9.5 


140 


G93; 800 


TY CrA 


19 01 40.79 


-36 52 34.2 


14-08-2005 


B9 


140" 


VOO: SOO 


T CrA 


i'l Oi ."-i.Tr! 


r,7 1').') 




FO,' 


i 10" 


FSl: soil 


VISIR N-band 


SY Cha 


10 56 30.45 


-77 11 39.3 


06-05-2006 


Ml) 


178 


K06; W97 


WX Cha 


11 09 58.74 


-77 37 08.9 


06-05-2006 


K7-M0 


178 


K06; W97 


HD 98922 


11 22 31.67 


-53 22 11.4 


07-05-2006 


B9 


>540 


A98; H78 


HD 101412 


11 39 44.46 


-60 10 27.9 


07-05-2006 


B9.5 


160 


T94; A05 


T Cha 


11 57 13.49 


-79 21 31.4 


04-05-2006 


G8 


66 


A93; A98 


WL 16 


16 27 02.34 


-24 37 27.2 


06-05-2006 


B8-A7 


125 


L99; G89 


SR 21A 


16 27 10.28 


-24 19 12.7 


06-05-2006 


G2.5 


125 


P03; G89 


Oph IRS 48 


16 27 37.19 


-24 30 35.0 


12-06-2005 


MO 


125 


O07; G89 


VV Ser 


18 28 47.86 


-1-00 08 39.8 


07-05-2006 


AOV 


269 


MOl; S96 


EC 82 


18 29 56.80 


-1-01 14 46.0 


07-05-2006 


MO 


269 


K06; S96 


NACO L-band 


WL 16 set 1 


16 27 02.34 


-24 37 27.2 


26-03-2006 


B8-A7 


126 


L99; G89 


WL 16 set 2 


16 27 02.34 


-24 37 27.2 


26-03-2006 


B8-A7 


126 


L99; G89 


WL 16 set 3 


16 27 02.34 


-24 37 27.2 


10-04-2006 


B8-A7 


126 


L99; G89 


Oph IRS 48 


16 27 37.19 


-24 30 35.0 


10-04-2006 


MO 


126 


O07; G89 



: all CrA sources assumed to be at same distance as that derived by SOO for HD 176386. 
References for spectral type, distance: A93: Alcala et al. (1993), A98: van den Anckcr et al. (1998), A04: Acke ^ van den Ancker (2004), AU04: Augereau 
& Papaloizou (2004), A05: Acke et al. (2005), C79: Cohen & Kuhi (1979), C06: Comeron (2007), D97: Dunkin et al. (1997), P84: Finkcnzeller & Mundt 
(1984), G89: de Gous et al. (1989), G93: Grady et al. (1993), G07: Grankin et al. (2007), H78: Houk (1978), H93: Hughes et al. (1993), H94: Hughes et al. 
(1994), KOI: Knude & Nielsen (2001), K06: Kcssler-Silacci et al. (2006), L99: Luhman & Rieke (1999), MOl: Mora et al. (2001), O07: Oliveira et al. (2007, 
in prep.), P03: Prato et al. (2003), R03: Rcsslor & Barsony (2003), SOO: Siebenmorgen et al. (2000), S96: Straizys et al. (1996), T94: The et al. (1994), 
VOO: Valenti et al. (2000), W97: Whittet et al. (1997), 



for the case of HD 100546 were the conditions good enough 
to match the source and standard spectra by applying dif- 
ferential airmass and seeing corrections. 

For the VISIR and NACO data, the extent of the spa- 
tial profile was derived as the standard deviation of the 
flux distribution Fi over the spatial pixels Xi. For each 
wavelength bin, the standard deviation cr^ was calculated 
as T,{xi-CYxFi/T,Fi, with centroid C = T,{xiXFi) /Y,Fi. 
Assuming a normal distribution, 3 a corresponded to 99% 
of the spatial extent of the flux. The NAOS-CONICA ob- 
servations of WL 16 and IRS 48 were taken using the adap- 
tive optics (AO) system. This resulted in high angular- 
resolution observations, which were typically close to the 
diffraction limit 0.1" at 3.3 fim). 



3. Results and discussion 

3.1. PAH detections and statistics 

The L-band spectra are presented in Figs. 1 and 2. Table 
2 lists the line intensities of the detected PAH features. 
The N-band spectra are shown in Figs. 3 and 4. 



We detect the 3.3 /j-m PAH feature in 6 out of 
24 sources, of which 4 Herbig Ae/Be stars (TY CrA, 
HD 100546, HD 98922, HD 101412) and 2 T Tauri stars 
(SR 21A, IRS 48). For VV Ser and T Cha, a marginal 
feature is found. The HI Br a, Pf 7, Pf S, and 10-5 transi- 
tion lines are seen in the majority of the sample, of which 
the HI Pf d line lies superposed on the 3.3 /xm feature. 
Note that we do not detect the 3.3 /zm feature toward the 
Herbig Ae star HD 141569. A tentative 11.2 fim feature 
was presented by Sylvester et al. (1996) and later clearly 
detected by Sloan et al. (2005). This absence of the 3.3 
and presence of the 11.2 /im feature is consistent with the 
model prediction of Li & Lunine (2003), where the 3.3 ijxa 
feature is weak due to the high degree of ionization in this 
disk (AUamandola et al. 1999). 

In the N-band spectra, PAH features are detected in 
3 Herbig Ae/Be (WL 16, HD 98922, HD 101412) and 4 
T Tauri stars (T Cha, SR 21A, IRS 48, EC 82), shown in 
Figs. 3 and 4. In all cases, except EC 82, both the 8.6 and 
the 11.2 ^m features are detected. For WL 16 and IRS 48, 
additional features at 11.9 and 12.8 /im are clearly seen. 

The 3.3 /Ltm PAH band is observed in only a very small 
number of T Tauri stars, 2 out of 18. Compared to Geers 



Geers et al.: Spatially extended PAHs in T Tauri disks 



5 



12 



10 




73 



o 4 



' I ' ' ' I ' ' ' I 

HI 10-5 Pf 6 
I I 



PAH 



I ' ' ' I ' ' ' I ' ' 

Pf 7 Br a 

I I 





3.0 3.2 3.4 3.6 3.8 4.0 
Wavelength (ytim) 

Fig. 1. ISAAC L-band spectra of sample with 3.3 /^m PAH 
detection. 



et al. (2006), all Herbig Ae and T Tauri stars with 11.2 /im 
Spitzer detections also show 3.3 /xm features in this L-band 
survey, confirming the detections. In addition, we find 3.3 
/xm PAH features in a few sources that were not included 
the Spitzer sample: IRS 48, WL 16, and HD 100546, all of 
which were known to have PAHs based on previous data 
(Malfait et al. 1998; Ressler & Barsony 2003; Geers et al. 
2005). Conversely, we have not found any new T Tauri 
disks with PAHs at 3.3 /im. The sensitivity varied sig- 
nificantly throughout the nights, but typical upper limits 
between 1 xlO-^^ and 5 x lO"!"^ W m ^ are obtained for 
the 3.3 /xm PAH feature, comparable to the Spitzer study 
of Geers et al. (2006) for the 11.2 /xm feature. The 11.2 
/xm feature is equally observed in only a small fraction 
of our T Tauri sample observed with VISIR. The Spitzer 
detections of Geers et al. (2006) are confirmed, and de- 
rived line intensities are consistent for most sources. The 
VISIR spectrum of VV Ser suffers from poor telluric cor- 
rection and does not show a clear 11.2 /xm feature, in con- 
trast to the Spitzer data. The spectra of HD 98922 and 
HD 101412 confirm the same weak, broad, and confused 
11.2 /im line noted to be possibly due to crystalline sili- 
cates (Geers et al. 2006). For both sources, Kessler-Silacci 
et al. (2006) confirmed the presence of crystalline silicates 
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from detection of the 23 and 33 /im crystalline silicate 
features. However, the 8.6 /im detections in both VISIR 
spectra strengthens the conclusion that PAHs also con- 
tribute to the 11.2 /im feature. In five sources, both the 
3.3 and 11.2 /xm PAH features are detected (HD 98922, 
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Fig. 4. VISIR N-band spectra of the sources with marginal 
(top 4) and no PAH detections (bottom 2). 



HD 101412, IRS 48, WL 16, SR 21A). One source, T Cha, 
has 11.2 /im clearly detected but 3.3 fim only marginally. 

3.2. Spatial extent 

Table 3 lists the radial spatial extent of the PAH fea- 
tures in half-width at half-maximum (HWHM). Spatial 
extent profiles are shown for a few sources, HD 100546, 
TY CrA, SR 21A, WL 16, and IRS 48, in the lower panels 
of Figs. 5-14. These figures show the spectrum in the top 
panel and in the bottom panel the diameter (full-width 
at half-maximum, or FWHM) of the spatial profile fit as 
a function of wavelength, for both the science source and 
the telluric and PSF standard star. 

The majority of sources with PAH detections show no 
spatial extent of the features beyond the surrounding con- 
tinuum emission from the disk at 3.3 /im, confining the 
source to the same spatial extent as the disk continuum 
emission. In many cases the upper limit on the spatial ex- 
tent is close to that of SR 21A, with a radius of 0.33" 
(41 AU at d = 125 pc), an example of a source with very 
strong unresolved PAH features (Fig. 11). These obser- 
vations constrain the origin of the PAH emission from T 
Tauri sources to the circumstellar disk. 

The hydrogen emission lines, 10-5, Pf S, Pf 7, and Br 
a, are detected in the majority of the spectra and are 
spatially unresolved in all cases. These lines are expected 
to originate from hot ionized gas in the very inner parts 
of the star-l-disk system, close to the surface of the central 
star. 

The maximum spatial extent can be observed either 
when the disk is fully face-on or when it is inclined at a 
moderate degree with the slit of the spectrograph aligned 
with the semi-major axis of the disk. For most sources, the 
inclination and position angle of the disk are unknown and 
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Fig. 5. (Upper panel) ISAAC L-band spectrum of 
HD 100546, scaled by a factor 1.25 to match the ISO spectrum 
and then shifted by +4 Jy for clarity. The grey band highlights 
the position and FWHM of the 3.3 PAH feature as deter- 
mined in this paper. (Lower panel) FWHM of the spatial pro- 
file, as extracted from 2D spectral image, of HD 100546 {black 
line) and its corresponding standard HR 4757 {grey line). The 
FWHM of the latter agrees with that of the HD 100546 con- 
tinuum after correction for airmass and seeing differences. 



thus the derived spatial extents are in most cases treated 
as a lower limit for the disk component. For a typical 45° 
inclination, assuming a circular disk and aligning perpen- 
dicular to the apparent semi-major axis, this would lead 
to an underestimate of the extent by a factor of 1.4. In 
the following, individual cases are discussed. 

3.2.1. HD 100546 

In HD 100546 (Fig. 5) the long-slit of ISAAC was ahgned 
to the semi-major axis of the disk, to measure the max- 
imum spatial extent of the emission. The 3.3 /im PAH 
feature is spatially resolved, with a radial extent of 0.40", 
found to be slightly larger than the neighboring 3 /im con- 
tinuum extent. After correcting for the difference in air- 
mass and seeing with respect to the standard star spec- 
trum, we find that the feature has a spatial extent of 
0.12" ± 0.032" and that the continuum emission is es- 
sentially unresolved. Assuming a distance of 103 pc, this 
corresponds to a radial extent of 12 ± 3 AU. 

HD 100546 is one of the most nearby isolated Herbig 
Be stars for which a circumstellar disk has been firmly 
established through space and ground-based spectroscopy 
and imaging. Physical parameters for HD 100546 and its 
distance have been derived from Hipparcos observations 
(van den Ancker et al. 1998): it is a B9Ve star with a 
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Table 2. Summary of PAH feature intensity in W m~^ and presence of HI lines. 



Source 


H Unes 


3.3 pm 


8.6 fim 


11.2 /^m 


12.7 fim 


SX Cha 


Y 


< 1.9 X 10~^^ 


- 


- 


- 


SY Cha 


Y 


< 8.0 X 10"^'' 


- 


- 


- 


WX Cha 


Y 


< 2.7 X 10"^" 


N 


N 


- 


HD 98922 


Y 


4.4 X 10~^* 


Y 


Y 


- 


HD 101412 


Y 


2.8 X 10~^^ " 


Y 


Y 


- 


HD 100546 


Y 


Y 


~ 


~ 




T Cha 


Y 


< 2.5 X 10~^^ 








IKAb 12b6o-7o26 


Y 


N 








HT Lup 


Y 


<■ l.O X lU 








GW Lup 


Y 


5; 2.4 X lU 








SZ 73 


Y 


9.4 X 10 








GQ Lup 


Y 


< 5.4 X 10"^'' 


- 


- 


- 


HD 141569 


Y 


< 6.9 X 10^^'' 


- 


- 


- 


IM Lup 


Y 


< 1.9 X 10"^-'^ 


- 


- 


- 


RU Lup 


Y 


< 5.7 X lO^-'^ " 






~ 


DoAr 24E 


Y 


< 4.1 X 10~^^ 








WL 16 


Y 


< 9.6 X 10"^^ 


Y 


Y 


Y 


Em* SR. 21A 


Y 


3.5 X 10"^-'' 


Y 


Y 


T 


Oph IRS 48 


Y 


2.4 X 10"^-'' 


Y 


Y 


Y 


Em* SR 9 


Y 


< 9.7 X IQ-i'* " 








Haro 1-17 


Y 


< 4.5 X 10^1'^ 








V1121 Oph 


Y 


1.2 X lO"^'^ " 








Wa Oph 6 


Y 


2.3 X lO"^'^ " 








VV Sor 


Y 


< 9.8 X lO-^'' 




T 




CoKu Ser G6 


T 


N 








EC 82 






N 


T 




HD 176386 


N 


< 4.2 X 10-^** 








TY CrA 


Y 


4.5 X 10-^^ 








T CrA 


Y 


< 3.9 X 10-^^ 









- : not observed 

" Line intensity affected by HI Pf 5 line 



Table 3. Summary of the radial spatial extent of the PAH feature (half- width at half-maximum). 



Source 


3.3 /^m 


8.6 




11.2 


IJ.ni 


12.7 




PAH detected and spatially resolved 


HD 100546 


0.12" 


12 AU 














TY CrA 


0.39" 


55 AU 














WL 16 






0.43" 


54 AU 


0.49" 


61 AU 


0.46" 


58 AU 


IRS 48 


0.11" 


14 AU 


0.25" 


31 AU 


0.32" 


40 AU 


0.33" 


41 AU 


PAH detected but spatially unresolved 


HD 98922 


< 0.24" 


< 130 AU" 


< 0.17" 


< 92 AU 


< 0.15" 


< 81 AU 






HD 101412 


< 0.24" 


< 38 AU 


< 0.16" 


< 26 AU 


< 0.15" 


< 24 AU 






Em* SR 21A 


< 0.33" 


< 41 AU 


< 0.16" 


< 20 AU 


< 0.17" 


< 21 AU 


< 0.15" 


< 19 AU 


T Cha 


< 0.34" 


< 22 AU 






< 0.19" 


< 13 AU 






EC 82 






< 0.18" 


< 47 AU 


< 0.23" 


< 60 AU 







: uncertain distance to source 



mass of 2.4 Mq. The circumstellar disk has been resolved 
through imaging at visible (Grady et al. 2001), near- 
infrared (Pantin ct al. 2000; Augcrcau ct al. 2001), mid- 
infrared (Liu et al. 2003), and millimeter (Wilner et al. 
2003) wavelengths. Augereau et al. (2001) observed an 
elliptical structure, consistent with an extended inclined 
disk, based on HST/NICM0S2 coronagraphic imaging at 
1.6 yum and derived values for the disk inclination and 
position angle adopted here. 

Bouwman et al. (2003) modeled the spatial distri- 
bution and chemical composition of the dust around 
HD 100546 and find that the data required the presence 
of a small grain component at '~200K, inconsistent with 
a uniform flaring disk (DuUemond ct al. 2001). Instead, 
they propose that the disk of HD 100546 has a gap be- 
tween 1 and 10 AU and is 'puffed up' at a radius of 



about 10 AU. Based on these results, they conclude that 
the observations of the spectral energy distribution (SED) 
of HD 100546 arc consistent with a largely clcarcd-out in- 
ner region between the inner rim and 10 AU and that 
most of the disk material is located in the outer parts of 
the disk. This conclusion was subsequently confirmed by 
Acke & van den Ancker (2006b). Our derived spatial ex- 
tent of the 3.3 /um feature of 12 ± 3 AU would place the 
PAHs responsible for the 3.3 /im emission at the puffcd- 
up rim at the outer edge of the cleared out region where 
a relatively large fraction of UV radiation is intercepted. 

Habart et al. (2006) report the presence of a 5-10 AU 
radial gap in the 3.3 /^m PAH spatial distribution, and a 
clear extension to at least 50 AU radially, which is larger 
than the typical extent (30 AU) found in their sample. 
This gap is not seen in our ISAAC data because of our 
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Fig. 6. (Upper panel) ISAAC L-band spectrum of TY CrA. 
(Lower panel) FWHM of spatial profile of TY CrA and stan- 
dard star HR 5812. 
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3.3 lAia 
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Fig. 7. Observed spatial emission profiles for WL 16 (left), 
IRS 48 (middle), and TY CrA (right), for the PAH feature at 
indicated wavelength (solid line) and continuum emission at 
10.51, 3.56, and 3.35 /xm respectively (dashed line). 



lower spatial resolution. Our ISAAC HWHM of the spa- 
tial extent of ~ 12 AU is smaller than their full 50 AU 
extent, presumably from our lower sensitivity to weaker 
extended emission with ISAAC. Their PAH distribution 
is consistent with our finding that the PAH emission is 
dominated by the dust rim at the outer edge of the gap. 

3.2.2. TY CrA 

For TY CrA (Fig. 6), the 3.3 jim feature is found to be 
spatially resolved (0.78", corresponding to a radial ex- 
tent of 54 AU) compared with the underlying 3.3 /xm disk 
continuum (0.73", or 51 AU radially). The coordinates of 
this compact component match within 1" with the optical 
multiple star system and with the nearest bright 2MASS 
Ks-band source. TY CrA is a very strong X-ray source, 
surrounded by a strong reflection nebula. Unpublished 
Spitzer/IRAC observations show a large extended nebu- 
losity at 8 /im around this source (Allen, priv. comm.). 
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Fig. 8. (Upper panel) VISIR N-band spectrum of WL 16. 
(Lower panel) FWHM of spatial profile of WL 16 and stan- 
dard star HR 6879. 
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Fig. 9. VISIR 2D spectral image of SR 21A (left) and WL 16 
(right); 1 pixel = 0.02719". 



confirming that PAHs are present on even larger scales as 
shown by Siebenmorgen et al. (2000) at 11.28 /im. The 
spatial emission profile of the 3.3 /im feature is shown 
in Fig. 7, which confirms the presence of the weak ex- 
tended PAH emission. PAH features have been previously 
observed toward this star with ISO (Klaas et al. 2006; 
Acke & van den Ancker 2006a). NACO imaging revealed 
this system to be a possible quadruple system of low-mass 
M type stars (Chauvin et al. 2003) , possibly affecting the 
circumstellar disk(s) structure and resulting PAH emis- 
sion. 

3.2.3. WL 16 

For WL16, the disk continuum emission at 8-13 /im is re- 
solved with a radial extent of about 0.4" (50 AU) (Figs. 8 
and 9). The 8.6, 11.2, and 12.7 /im PAH features are re- 
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Fig. 10. (Upper-left panel) NACO L-band spectrum of 
WL 16, slit parallel to semi-major axis of the disk. (Lower- 
left panel) FWHM of spatial profile of the "parallel" obser- 
vation of WL 16 (black) and standard star HD 142378 (grey). 
(Upper-right panel) NACO L-band spectrum of WL 16, slit 
perpendicular to semi-major axis of the disk. (Lower panel) 
FWHM of spatial profile of the "perpendicular" observation of 
WL 16 (black) and standard star HD 142378 (grey). 



solved with respect to the continuum, with radial spatial 
extents of 0.43" (54 AU), 0.49" (61 AU), and 0.46" (58 
AU) . The spatial emission profile is shown in the left panel 
of Fig. 7. The PAH emission at 10-13 fim has broad wings 
in the spatial direction, extending beyond the spatial ex- 
tent derived using the statistical derivation described in 
Sec. 2.5, illustrating how in certain cases it is hard to 
quantify the spatial extent in a single number. The extent 
varies with wavelength, which is not expected from uni- 
form diffuse background PAH emission. Deriving the ex- 
tent alternatively as 1% of the peak level, we find a value 
of 80-140 AU, which is about a factor 2.5 lower than that 
found by Ressler & Barsony (2003) of 440 by 220 AU ra- 
dially, based on Keck imaging covering 7.9 to 24.5 fim. 
Given the reported semi-major over semi-minor axis ratio 
of 0.466, the effect of non-alignment of the VISIR slit with 
the semi-major axis can be at most a factor of 2. 

NACO L-band spectroscopy was taken in two settings, 
aligning the long-slit parallel to the semi-major axis and 
perpendicular. The 3.3 /xm feature is undetected in both 
orientations (Fig. 10). The radial spatial extent of the con- 
tinuum at 3.3 fim is small and very similar between both 
orientations. The non-detection at 3.3 /xm would be con- 
sistent with predominantly ionized PAHs (Allamandola 
et al. 1999; Li & Draine 2001). 

3.2.4. SR 21A 

In contrast with WL 16, SR 21 A is a source with very 
strong 3.3, 8.6, and 11.2 fim features but spatially unre- 
solved. The 3.3, 8.6-11.2 /xm disk continuum emission has 
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Fig. 11. (Upper panel) ISAAC L-band spectrum of SR 21A. 
(Lower panel) FWHM of spatial profile of SR 21A and stan- 
dard star HR 6629. 



a radial extent of < 0.15-0.33", corresponding to ~19-41 

AU, shown in Figs. 11 and 12, and in the left panel of 
Fig. 9. This source has been interpreted as a 'cold disk' 
source by Brown et al. (2007) with a gap between 0.45 
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Fig. 12. (Upper panel) VISIR N-band spectrum of SR 21A. 
(Lower panel) FWHM of spatial profile of SR 21A and stan- 
dard star HR 6879. 
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and 18 AU, but with a large outer disk. The origin of the 
PAH emission is suspected to arise at the outer edge of 
the inner disk or within the gas filling the gap, inside 30 
AU radius. The unknown position angle may have led to 
a chance alignment of the slit close to the semi-minor axis 
of the disk, which may have reduced the apparent spatial 
extent. 

3.2.5. IRS 48 

IRS 48 is the only source in this sample for which all PAH 

features are spatially resolved. IRS 48 is a low-mass MO 
star in the p Oph cloud for which strong 3.3, 8.6, 11.2, 
11.9, and 12.7 /itm PAH features are detected. For this 
source VISIR images have shown the presence of a 30 AU 
radius gap in the inner disk (Geers et al. 2007) with the 
PAH feature emission originating from a region inside the 
gap. The disk continuum emission at 3-4 /im in the NACO 
spectrum is spatially resolved with a radial extent rising 
from 0.06-0.09" (7.5-11 AU) with wavelength. The 3.3 
/im PAH feature is detected and spatially resolved on a 
radial scale of 0.11" (or 15 AU) above the extent of the 
continuum of 0.07" or 9 AU (Figs. 7 and 14). Assuming 
diffraction limit and no influence of seeing and correcting 
for the contribution of the PSF based on the spatial extent 
of the underlying continuum, the radial extent of the 3.3 
/xm features is found to be ^ 11 AU. The disk continuum 
emission at 8-13 /xm in the VISIR spectrum is resolved 
with a radial extent increasing from 0.25"-0.33" (31-41 



AU). The 8.6, 11.2, and 12.7 /zm features are spatially re- 
solved on a radial scale smaller than the surrounding disk 
continuum, with radial extents of 0.25, 0.32, and 0.33" 
(31, 40, 41 AU) respectively, half-width at half-maximum 
(HWHM). These extents arc consistent with the full (99%) 
radial 8.6 and 11.2 pm PAH extent of 75 and 90 AU re- 
spectively, found in spatially resolved in VISIR images in 
Goers ct al. (2007). Also, our finding of the smaller feature 
extent compared with the continuum is consistent with the 
VISIR images showing PAH emission to be smaller than 
the mid-IR 18.7 fim emission. 

3.2.6. T Cha 

For T Cha, van den Ancker et al. (1998) measured a dis- 
tance of only 66 pc, making it one of the closest, isolated 
T Tauri stars, and not part of the Chamaeleon cloud 
[d = 178 pc). A weak 11.2 /xm feature and a marginal 
3.3 jim feature are detected. The 3.3 jim feature is un- 
resolved, with a radial spatial extent of the continuum 
emission of <0.34", corresponding to < 22 AU. The 11.2 
fxm. feature is also unresolved, with a continuum spatial 
extent of <0.19", corresponding to a radial extent of <13 
AU. 



3.2.7. Summary and discussion 

Table 3 summarizes our detections of spatially re- 
solved PAH features. It tabulates the half- width at half- 
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Fig. 15. Radial extent of 3.3 fj,m feature versus 11.2 fim fea- 
ture; Herbig Ae stars indicated by filled diamonds, T Tauri 
stars by filled triangles; models with A'c = 96 indicated by 
grey, A^c=50 by black. Observations of IRS 48, SR 21 A, T Cha, 
HD 101412 and HD 98922 are indicated by open symbols. 
Model prediction point extracted from Fig. 4 of Habart et al. 
(2004) for a disk around a Herbig star for A^c = 100 is indicated 
by a cross. 

maximum extent of the emission, assuming it can be fit- 
ted with a Gaussian profile. Under this assumption, this 
HWHM extent corresponds to the radius of the disk, in- 
wards of which 76% of the normalized cumulative emission 
originates. 

This HWHM extent is compared with model predic- 
tions by Visser et al. (2007), reproduced in Fig. A.l 
in Appendix A, with an improved extraction technique. 
Figure 15 plots the radial extent of the 3.3 /im versus 
the 11.2 /xm feature, for both model predictions of Herbig 
Ae stars and T Tauri stars, for both small (7Vc=50) and 
large (iVc=96) PAHs. The model extent of the 3.3 /im fea- 
ture varies strongly with PAH size, decreasing by a factor 
2-3 when the PAH size is increased by a factor of ^ 2, 
while the extent of the 11.2 /im feature is less sensitive. 
For Herbig Ae stars, the model predicts for small PAHs 
(A^c=50) that ^76% of the PAH emission should originate 
from a ring with a typical radius of ^160 AU for 3.3 and 
11.2 /im and 175-190 AU for 6.2, 7.7, and 8.6 fim. For 
larger PAHs (A^c=96) the predicted distribution is more 
centrally peaked, with 76% of the emission arising inwards 
of radii of 47, 96, 105, 103, and 117 AU for the 3.3, 6.2, 7.7, 
8.6, and 11.2 /tm features respectively. In T Tauri disks the 
small PAH emission is also more centrally peaked, within 
radii of 15 AU for the 3.3 /tm feature, and 50-60 AU for 
the 6.2, 7.7, 8.6, and 11.2 /xm features. 

The observed extents for IRS 48, SR 21A, T Cha, 
HD 101412, and HD 98922 are included in Fig. 15. Most 
PAH detections have upper limits on the radial spa- 



tial extent of typically 20-40 AU (with the exception of 

HD 98922, whose distance is uncertain) while for a few 
sources the PAH emission is spatially resolved with typi- 
cal extent of 12-55 AU. One source, IRS 48, has spatially 
resolved detections of all features. 

For all observations with detected PAH features, the 
(upper limit for the) radial extent of the 11.2 /im feature 
is a factor of 1.5 6 lower than derived from the prediction 
of the template model and we do not see any evidence for 
ring-shaped PAH emission as the model predicts. Within 
the Visser et al. (2007) models, this implies larger PAHs 
of typically 100 or more carbon atoms. The model predic- 
tion by Habart et al. (2004) for Herbig Ae stars with large 
PAHs (A^c=100) is also indicated in Fig. 15. Our measure- 
ment for the Herbig Ae star, HD 101412, is a factor 6-7 
smaller than their and our prediction. 

The typical extent of the 3.3 /im feature at 76% of 12- 
55 AU agrees with the predictions of both Visser et al. 
(2007) and Habart et al. (2004). Our derived upper lim- 
its on the extent of the 8.6 and 11.2 /im features of typ- 
ically < 30 AU, are about a factor 2 smaller than the 
model predictions, which also hints at the presence of 
larger PAHs than considered in those models. As noted 
by Geers et al. (2006) , a significant number of these PAH 
sources are known to have substantial gaps out to 40 AU 
in their disks (Brown et al. 2007). The presence of warm 
3 /tm disk emission suggests these sources have gaps, not 
holes. These observations put the PAH emission in the in- 
ner part of the disk, and in some cases (i.p. HD 100546) 
it very likely comes from the dust at the outer edge of the 
gap. 

Four out of the five sources with spatially unresolved 
PAH features (leaving out HD 98922) have average upper 
limits on the extent of 30-40 AU, with no apparent dif- 
ference between the Herbig Ac star HD 101412 and the 
3 T Tauri stars. Potential misalignment of the slit with 
the semi- major axis of an inclined disk will reduce the 
measurable spatial extent (see also Appendix A) . In addi- 
tion, unresolved features could have at least two geomet- 
rical causes. A larger inclination of the disk (from face- 
on) would decrease the measurable extent, although at 
near-IR wavelengths, for close to edge-on disks, scattering 
could dominate the spatial extent of the source of emis- 
sion. Four out of the five unresolved sources (excluding 
EC 82) have no indications of a strong inclination in their 
SED (Kessler-Silacci et al. 2006). The second geometrical 
effect is the flaring or the covering fraction of the disk. 
For a disk with a decreasing covering fraction (through 
decreasing scale height), the feature to continuum ratio of 
the PAH features will decrease for very flat disks (Geers 
et al. 2006). Both HD 101412 and HD 98922 have SEDs 
consistent with a flat type II source (Meeus et al. 2001), 
while the SR 21A and T Cha both have SEDs consistent 
with still moderately large scale height, although the SED 
interpretation is more complicated due to the suspected 
presence of a gap (Brown et al. 2007). Besides geometri- 
cal and instrumental effects, the alternative physical in- 
terpretation for the small extent could be the presence of 



12 



Geers et al.: Spatially extended PAHs in T Tauri disks 



predominantly large PAHs, which model predictions show 
to radiate mostly from smaller scales, as discussed above 
(Visser et al. 2007). 

4. Conclusions 

We detect the 3.3, 8.6, and 11.2 /xm PAH features in a 

small fraction of our sample of T Tauri stars, with typical 
upper limits between 1 xlO~^^ and 5 x 10"^'' W m~^. 
Compared with the Spitzer survey in Geers et al. (2006), 
wc confirm all their 11.2 /im detections and find marginal- 
to-clear 3.3 /um detections. The bias in our sample pre- 
vents us from deriving an independent detection fraction, 
but the low number of PAH observations toward T Tauri 
stars is consistent with the low detection rate of 8% by 
Geers et al. (2006). In two specific sources, WL 16 and 
HD 141569, the PAHs are expected to be largely ionized 
because of the strong 7.7 and 8.6 /xm features and lack of 
3.3 /zm. 

The spatial extent of the PAH features is shown to be 
confined to scales smaller than 0.1 — 0.4" (HWHM), cor- 
responding to radial scales of 12-60 AU in the disk, at 
typical distances of 150 pc, barring exceptional cases with 
unknown or much larger distances. In a few examples, the 
PAH features are resolved to be more extended than the 
hot underlying dust continuum of the disk, whereas in one 
case, IRS 48, the extent of the PAH emission is confirmed 
to be less than that of the large grains. The typical ex- 
tent of the PAH features of 15-60 AU is found to be very 
similar for both Herbig Ae and T Tauri stars, and simi- 
lar for all features. For Herbig Ae stars, the small 12-55 
AU extent and absence of any ring emission is consistent 
with the model predictions of larger (>100 carbon atoms) 
grains. The same conclusion of a need for large PAHs holds 
for the T Tauri stars where the 8.6 and 11.2 /im features 
appear to be smaller than predicted, although here the 3.3 
/xm extent is consistent with smaller (^ 50 carbon atoms) 
PAHs. Given the large fraction of disks with gaps and 
PAHs, future modeUng studies of the PAH extent should 
include the presence of gaps in disks. 
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Appendix A: Spatial extent models 

Visser et al. (2007, hereafter V07) modeled the chemistry 
of and infrared emission from PAHs in circumstellar disks, 
including an analysis of the spatial extent of the emis- 
sion. The algorithm to determine the spatial extent of the 
emission was further improved to allow for a more direct 
comparison with observed spatial profiles by summing the 
emission in a (narrow) slit rather than in concentric rings 
as done in V07. A simplification was introduced by mea- 
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Fig. A.l. Normalized cumulative intensity of the five main 
PAH features (black) and the continua at 3.1 and 19.6 /xm 
(gray) for CsoHig and C96H24 in the standard Herbig Ae/Be 
and T Tauri model disks (including PAH chemistry) from 
Visser ct al. (2007). The disks were synthetically observed facc- 
on through an infinitely wide slit. The vertical bars indicate 
where each curve reaches 76%, corresponding to the HWHM 
of the Gaussians in Sect. 2.5. 



suring continuum-subtracted peak intensities instead of 
spectrally integrated feature fluxes. Since the model fea- 
tures all have the same shape (Draine & Li 2007), this 
produces the same spatial profiles. The new procedure is 
essentially the same as outlined in Sect. 2.5 for the ob- 
servations. The improved profiles for the standard Herbig 
Ae/Be and T Tauri models (i?disk,in = 0.077 AU, i?disk,out 
= 300 AU, Mdisk = 0.01 M0) from V07, synthetically ob- 
served face-on through a slit wide enough to cover the 
entire disk, are presented in Fig. A.l. 

If the source docs not fit inside the width of the slit, 
part of the emission originating at large radii is blocked, 
thus generally reducing the apparent spatial extent com- 
pared to observations through a wider slit. This is shown 
in Fig. A.2. With a 40 AU slit (0.27" at 150 pc), the loca- 
tion inside of which 76% of the feature emission originates 
shifts to smaller radii by the order of 10 — 15% compared 
to an infinitely wide slit. 

The effects of inclination are less straightforward. For 
objects at an inclination of about 45°, the observed spatial 
profiles differ only slightly from the face-on profiles if an 
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Fig. A. 2. Normalized cumulative intensity for CsgHig in the 
standard Herbig Ae/Be model disk (including PAH chemistry) 
from Visser et al. (2007). The disk was synthetically observed 
at 45° or through a narrow slit as indicated. Lines and colors 
are as in Fig. A.l. 

infinitely wide slit is used. With a narrow slit, part of the 
emission is obscured asymmetrieally across the width of 
the slit and the emission can appear to become both more 
or less extended, depending on the details of the source. 
This also occurs when the source is viewed more edge on, 
regardless of the width of the slit. 



